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Abstract 

In order to investigate more features of the Brans-Dicke cosmology in 
the five-dimensional space-time, we explore the solutions of its dynam- 
ical systems. A behavior of the universe in its early and late time by 
means of the scale factor is considered. As a results, we show that it is 
possible to avoid the big rip singularity and to cross the phantom divide 
line. Furthermore, we review the dark energy component of the universe 
and its agreement with the observation data for this 3-brane Brans Dicke 
cosmology by means of the cosmological parameters. 

1 Introduction 

The Brans-Dicke (BD) theory [1], the simplest case of the scalar tensor theory 
[2], is defined by a constant coupling parameter u> and a scalar field 4>. General 
relativity is recovered when u goes to infinity [3] and from timing experiments 
using the Viking space probe [4], uo must exceed 500. This constraint ruled out 
many of extended inflation [5, 6] and provides a succession of improved models 
of extended inflation [7, 8, 9, 10]. Furthermore, all important features of the 
evolution of the Universe such as: inflation [11], early and late time behavior of 
universe [12], cosmic acceleration and structure formation [13], quintessence and 
coincidence problem [14], self -interacting potential and cosmic acceleration [15], 
High energy description of dark energy in an approximate 3-brane [16] could 
be explained successfully in the BD formalism. For a large value of the u- 
parameter, BD theory gives the correct amount of inflation and early and late 
time behaviors, while small and negative values explain cosmic acceleration, 
structure formation and coincidence problem. 

The dark energy, qualified as responsible for the cosmic acceleration deter- 
mines the feature of a future evolution of the universe. The nature of this kind 
of energy may lead to the improvement of our pictures on particle physics and 
gravitation. The investigations on the nature of the dark energy lead to various 
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candidates. Among them, the most popular one, the cosmological constant [17], 
the dynamical scalar fields, described by the equation of state p = (7 — l)p, 
like quintessence [15, 18, 19, 20] (7 > 0) or like phantom [21, 22, 23] (7 < 0). 
However, many problems associated with the phantom candidate have to be 
explained. The first one is the big rip [24], the appearance of a future sin- 
gularity in a finite time, at which all cosmological parameters blow up. To 
overcome this problem, many models, whose solutions do not suffer from big 
rip, have been proposed such as the model which stated that the final state 
of phantom cosmology may be inflation [25], the one inspired from the string 
theory [26], and the model in which the de-Sitter solution is attractor of phan- 
tom cosmology [27]. The second one is how to cross the phantom divide line 
7 = 0. However, the dynamical transition from 7 > to 7 < or vice versa 
is possible [28, 29, 30, 31]. We are thus motivated to explore the cosmological 
behavior of the phantom field, in such a way that inflation is the generic feature 
of phantom cosmology rather than the big rip as in Ref. [25] and also to explain 
the possible crossing of phantom divide line. 

More recently, a great deal of interest is devoted to high dimension space- 
time. Superstring, which predicts a new type of nonlinear structure, called a 
brane [32, 33, 34], suggests that our universe might be of higher dimensions 
which are compactified. The matter field content of our universe is confined to 
a four dimensional space-time called a 3-brane in the case of (5D) space-time. 
In addition, several works have studied higher dimensional BD theory in order 
to use the advantages of the combination of the high dimensional cosmology 
and the BD theory. 

In a recent paper [16], we have generalized the BD cosmology to a 3-brane 
with a nonzero cosmological constant, A4, derived from 5-dimensional (5D) 
bulk space-time. In this paper, we discuss the evolution of the scale factor 
at the very early stage of the expansion of the universe and at the late stage, 
where a new kind of matter (dark energy) is suspected to be responsible for 
the cosmic acceleration. We show that a possible accelerated era could appear 
if the universe undergoes a bounce state in the past, and avoids the big bang 
singularity or a turnaround state avoiding consequently the big rip singularity. 
Otherwise, the accelerating era appears with this two singularities in the past 
or/and in the future. Furthermore, in order to generalize our previous work [16] 
and to complete this study to late time accelerating universe, we linearize the 
dynamic system, by means of the cosmological parameters, in the intermediate 
energy levels. 

This paper is organized as follows. In section 2 we present the equations of 
the fields in the 3-brane Brans-Dicke cosmology. These expressions will be given 
in the shape of a dynamic system by introducing some variables. In section 3, 
we discuss the evolution of the scale factor and the possibility to overcome 
both the big rip singularity and the crossing divide line in phantom cosmology, 
by assuming that the universe started out small. In section 4 we examine 
the late time accelerating universe by giving the solution in the intermediate 
energy levels around the stable equilibrium solution. Section 5 is devoted to 
the conclusions. 
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2 Field Equations in 3-brane Brans-Dicke cosmology 



The detail of the 3-brane Brans-Dicke cosmology derived from 5-dimensional 
bulk space-time is given in [16]. In this section we shall review the main results. 
First, we consider that the behavior of BD field is sensitive only to a physical 
3-brane. Thus, it is described by the same action as in 4-dimension (4D) 

S = J d 4 x^9 (V [R - 2A 4 ] - - 16Tr£ m ^ . (1) 

Furthermore, to recover the BD cosmology at low energy, we add simply a 
BD stress energy tensor to the Einstein equation, G^ v [32], where all quadratic 
and mixed terms of this stress-energy tensor are canceled. Therefore, the mod- 
ified Einstein equation takes the form: 



A G^ U = -A 4 ^ + 8ttG n t^ + fc|n^ (2) 

+ ^(fci* fc>> -2<W;a (t> ]X ) + ^(<?W -Qnu^4>) 

where A4 is the 4D cosmological constant, 8itGn = k\ is the 4D gravitational 
constant (G/v is the Newton's constant of gravity), k\ is the 5D gravitational 
constant and t^, Yl^ v are respectively the energy momentum, the quadratic 
tensors on the brane. 

We have shown that the equation of motion of the BD field and the equation of 
state are modified. Indeed, by varying the action (1) versus the metric tensor 
and BD field (p, in the homogeneous and isotropic Friedman-Robertson- Walker 
metric with scale factor a(t) and spatial curvature index k (with Gn = ^), one 
obtain the following BD field equations: 



d A4 47F /q n\ k ht n 2 u(4>\ 2 la4> 10 



a 2 k A4 8ir k* 2 uj ( <p\ a0 
a 2 a 2 3 3</>^ 36^ 6 \<t>) a<f> 



p = -3^(p + p) (6) 

where p and p are respectively the energy density and the pressure of the cosmic 
fluid with the equation of state p = (7 — l)p. 

Equation (5) is the modified BD field equation, while in the standard BD cos- 
mology, it can be written as (uj 7^ — |), 

ld ^- 87T -((3 7 -4)p-^) (7) 



a 3 dt 3 + 2w vv ' Jr 4vr 
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Comparing the BD field equations (3)- (5), we conclude that at high energy 
limit (p 2 3> p), the 3-brane BD cosmology is described by the same manner 
as in the 4D BD cosmology with the equation of state of a perfect fluid is [16] 
p=(2 7 -l)p. 

By introducing the variables, H = |, F = | and Z = the field equa- 
tions (4), (5) and (6) become: 



dH 
~~dt 



-H 2 - |F 2 



+ HF 



2lo 



:A 4 



3(3 + 2uj)' 

(37W - 2uj + 3) (67W -2^ + 3) /cf 
"TT^ Z , <, \ TT7T; 



(3 + 2u) 



(3 + 2u) 36 ' 



(8) 



dF _ F 2 o ffF 1 4 (37 ~ 4) (3 7 -2) fc| 2 



(9) 



= H — —F + iTF H — ^ - % - T^P (10) 

6 cr 3 36 

3 Scale factor evolution in 3-brane Brans-Dike the- 
ory 

In this section we use the preceding dynamic equations to discuss the evolution 
of the scale factor in three cases: low, high and intermediate energy levels. We 
suppose that near the big-bang time, the universe is located in the equilibrium 
state, as in Ref. [35], namely a = 0. 



3.1 Low energy case: 

When the energy density p is sufficiently diluted in the universe, we are prac- 
tically in low energy limits, i.e. /;> p 2 . The preceding field equations (8), (9) 
and (10) become by eliminating the Z- variable: 



(2cj + 3)^ = -3H 2 (juj + 2) + ^F 2 uj(-fuj-2uj-l) + FHuj(4:-3-f) 

(XL £ 

+A 4 (70; + 1) - (3 7 w - 2co + 3) 4 (11) 



dF 1 
(2cj + 3)— = -3(3 7 - 4)# 2 + -F 2 (3~/lu - 8lu - 6) - 3FH (3 7 + 2lo - 1) 
ctt 2 

+A 4 (3 7 - 2) -3(37-4)A (12) 

cr 

which are exactly the field equations in the 4-dimensional BD cosmology. 
For a static universe (a = 0), the field equations (11) and (12) have a 
solution 
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F = a= 1™* (13) 
V (2w + 3) V ; 

and the BD field behavior depends on the 4D cosmological constant 

(t> = fae at (14) 

and, for 1 u> ^ — 1 



a = a* = J — 4-^ — tt- = — \ -, -r with k = l (15) 

Y A 4 (lo + 1) a]j (w + 1) 

Which means that in a static universe, the BD field evolves exponentially 
as e . 

Let us notice that the two solutions 4>(t) and a(t) are expected to be stable 
and we refer to [35] where the stability of such solutions is well established. 

Now, we investigate how the matter content of the universe, with an equa- 
tion of state p = (7 — l)p, affects its behavior compared to this stable solution. 

In this aim, we keep the BD field varying as e at and the previous equations 
become: 



dH 1 

(2lj + 3)— = -3H 2 (-fuj + 2) + -a 2 uj{joo-2uj-l) + aHuj{4-3-r) 

k 

+A 4 (7W + 1) - (3 7 w - 2uj + 3) — (16) 



= -3(3 7 - 4)H 2 + ^a 2 (3juj - 8lo - 6) - 3uH (3 7 + 2u - 1) 

+A 4 (3 7 -2)-3(3 7 -4)^. (17) 

Or, by eliminating the term 

3(37 - 4)— = -3H 2 (3 7 - 4) + 3aH (3 7 - 2w + 3 7 u; - 1) , (18) 
we obtain the following equation 

(37-4)^ = 0^(37-2^ + 370;- 1). (19) 

If the universe is dominated by the radiation (i.e., 7 = |), we recover the 
static universe and it remains eternally. Otherwise, we rewrite this equation, 
by putting = d, in the form 

e (37-4) y 1 



x the case uj = —1, which gives an eternally flat and static universe, is excluded in this 
study. 
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which gives easily 



6 = 9*e a BF=aj {t ~ U \ (21) 

where 9* is a constant of integration. Hence the scale factor varies as: 

o(t ) ^eW-M+c (22) 



where (3 = J i^t-w ^ 7 J^~ir ^ an< ^ c ^ s * ne integration constant. 



(2^+3) (3 7 -4) 

The behavior of the universe depends on the product (39*. Indeed, if the 
universe undergoes an era in which 9* is positive, it begins to accelerate for 
(3 > 0. The condition, in wich 9* is positive, is a feature of an era following the 
one wich we suppose to be a bounce state caracterized by a(t) — > and d(t) > 
[36]. Otherwise, the accelerating era is for (3 < 0, wich can be a feature of the 
era preceding a turnaround state caracterized by a(t) — ► and a(t) < [37]. 
The time in which 0* > could happen in the past, as in the high energy limit 
(see the next section) or in the future as in the low energy limit. In the former, 
the existence of the bounce state means that the big bang singularity could be 
avoided. In the future, however, it is the turnaround state which is important 
since it may avoid the big rip singularity if never it exists. 

To illustrate this situation, we consider the case uj » 1 i.e. (3 = y/Ki^^uo^ . 

This 



-oo, § 



U 



f , +cx> 



ike dust or radiation) and 



For (3 > 0, the parameter 7 varies in the range 7 € 
range excludes the ordinary forms of matter /energy ( 
the one like a cosmic string, 7 = |, in which the universe remains static and 
eternally. In other word, the accelerating universe caused by an exotic form 
of matter /energy is expected in the interval range where 7 < | (p < — |p) 
or 7 > |. In the case 7 < |, a possible nature of the dark energy will be a 
quintessence, (7 > 0), a domain wall, 7 = |, a cosmological constant, 7 = 0, 
or a phantom, (7 < 0). While for 7 > | the only possible candidate is the 
quintessece . In the case where (3 < 0, the accelerating universe is caused by 
the dark energy component of the energy density of the universe, for 7 £ §, § , 
and the quintessece is the only possible candidate for the nature of this kind of 
matter /energy. 

3.2 High energy case 

in the Ref. [16], We have shown that, at high energy limit, the 3-brane BD cos- 
mology could be described by the 4D BD cosmology with the following equation 
of state p = (27 — l)p, i.e. the 7-parameter in the standard cosmology is equal 
to twice the one in 5D bulk space-time. 

The field equations in this limit are the same as (11) and (12) where the 
7-parameter should be replaced by 27. 

The solutions, in this case, are the same as the one in the low limit with the 
equation of state p = (27 — l)p and (3 = ^ (2^+3) < ' 67 ~2 2 (37-2) J ~ 1 ' > 

We can draw the same conclusion by replacing the 7-parameter by 27. In- 
deed, for (3 > 0, the parameter 7 varies in the range 7 G 



-00, \ 



U 



§, +00 
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for to S> 1. This range excludes the forms of matter /energy like a domain wall, 
7 = ^, a cosmic string; 7 = |, or a quintessence, for which 7 6 |, | . In other 
word, the accelerating universe is caused by an exotic form of matter /energy, 
like a phantom or a cosmological constant, for 7 < | in agreement with [38] or 
like a quintessence for 7 > |. Otherwise, the cosmic acceleration is caused only 

by a dark energy, for which 7 £ 3, § , and the only possible candidate for the 
nature of this kind of matter/energy is, once more, the quintessece. 

3.3 The intermediate energy level 

In the intermediate energy level case, where both p and p 2 contribute to the 
evolution of the universe, the field equations (8), (9) and (10) become: 

(3 + 2co)^ = -6 (7^ + 1) H 2 - co ( - (2 ~ 2 27) + - ) F 2 - u;(6 7 - 4)ffF 
k 

-^(w(6 7 -2) +3) +A 4 (2 7 u; + 1) +7wZ (23) 



dF 

(3 + 2lj)— = -3(67-4)tf 2 -(w(4-37)+3)F 2 -3(67-l + 2w)#F 



dt 



3k 

2 (67 - 4) + A 4 (67 - 2) + 3 7 Z. (24) 



Since, in the low and high energy limits, the BD field increases as e at , we 
suggest that, in the intermediate energy limit, this BD field behaves in the same 
manner i.e. as e at and letting, only, the scale factor to be variable. With this 
assumption the equations (23) and (24) become: 

(3 + 2u,)^ = -6 ( 7 co + 1) H 2 - co ( - (2 - 2l) + 1 ) a 2 - ^(67 - 4)ff a 
- 4 (w (67 - 2) + 3) + A 4 (2 7 w + 1) + -jcoZ (25) 



= -3(67-4)# 2 -(w(4-37)+3)a 2 -3(67-l + 2w)£fa 
3A; 

- ^2 (67 - 4) + A 4 (67 - 2) + 3 7 Z (26) 
or by eliminating the Z-term between the two latter equations, one finds: 

3 (da + d 2 ^ = 3coada — 3k + ^-cua 2 + A 4 ^ a 2 (27) 
which becomes in term of a new variable 9 = a 2 : 

e-coaO- 2A 4 CJ + 1 - 9 = -2. (28) 
2co + 3 

We notice that this dynamical equation, and therefore its solution, is free of 
7— parameter and all kinds of matter/energy are welcomed. Its solution is: 
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9 = a 2 = al + c ie ~ at + c 2 e Q ^ +1 )' (29) 
where a* = ^^TJ for w / -1,. 

One can notice that at i — > the scale factor varies as ~ \ft and at t > but 
not too much, the evolution of the scale factor, described by a(t) 
consistent with the primordial rapid inflation. 

we conclude that, the 3-brane BD cosmology at the intermediate limit where 
p and p 2 are both considered, we recover the standard like expansion of the 
universe for all kind of matter /energy, at t — > 0, as in the standard cosmology 
but for radiation era only. For a late time t > 0, we recover an exponential 
expansion for all kind of matter /energy, however it is only for vacuum energy 
in the standard cosmology. 

We stress that the last two appriximations, denoted by low and high energy 
limit, should be a good and simple illustration to overcome the big rip problem 
and may explain how to cross the divide line 7 = in phantom cosmology. 
Indeed, for the big rip problem, if we replace the equation of state parameter 
7 by 7 — I in the low energy limit, i.e. we exclude the radiation from the 
matter/energy content of the present universe which is well justified by obser- 
vation [39], we conclude the following: The low energy limit corresponds to 
the late time evolution of the universe and the high energy limit to the early 
time. For the late time, we assume that the acceleration occurs in the case 
when f3 < while for the early time we assume that it happens in the case 
(3 > 0. Hence the early time acceleration is caused by the matter/energy con- 
tent with 7 G —00, 1 U ]2, +oo[ while for the late time the acceleration is due 

to the matter /energy content with 76 — 1,0 . We notice the possibility of a 
transition from a state of the universe, dominated by either a dust (7 = 1), a 
cosmological constant (7 = 0), or a quintessence (7 G]0, |[U]2, +00 [) to a state 
of the universe, dominated by a phantom (7 < 0), but undergoing a turnaround 
state avoiding therefore the big rip singularity. Hence, the final state of phan- 
tom cosmology may be inflation rather than big rip since the turnaround state 
stops the acceleration. 

Furthermore, it is also possible to cross the phantome divide line (7 = 0). 
Indeed, if the universe, at early time, is described by a matter/energy content 
as dust, cosmological constant, or quintessence (7 € 0, | ), the acceleration 
of the universe at late time is caused necessarly by a phantom like field whith 



7 G 



■2 



i.e. a transition from 7 > to 7 < 0. 



4 Late time accelerating universe 

In this section, we consider the intermediate case in which p and p 2 are both 
considered. The limit case, low and high energy, were considered in our previous 
work [16]. We have shown that our results are in agreement with the observation 
data, more precisely with the dark energy via the cosmological parameters. 
We can analyze how much today's the universe is far from late-time inflation 
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by linearizing the dynamical system about the stable cosmological solution with 
flat space and show how the Hubble parameter varies with the scale factor a(t). 



4.1 Stability 

Since the combined results of the cosmic microwave background and type la 
SNe [40, 41, 42] conclude that the universe undergoes a flat period today, we 
neglect the curvature parameter k/a 2 as a(t) increases with the expansion of 
the universe. Under these considerations, and in analogy with the previous 
section, the stable solution for (23) and (24) is: 



/ 2A 

' g " F "^ = V (^ + 3)(k, + 4) '" + 1 - 1 -'"- (30) 

Indeed, the stable solutions are obtained from the equilibriums points. To 
this end, we add to (23) and (24), the equation 2 

^ = -(^H + F)Z. (31) 

Neglecting the fc-term, the equations (23) and (24) at the equilibrium points 
become 



= -6 ( 7 u; + 1) H 2 - to ( - (2 2 27) + - ) F 2 - u;(6 7 - A)HF 

+ A 4 + 1) + jojZ (32) 

= -3 (67 - 4) H 2 - (uj(4 - 37) + 3) F 2 - 3(67 - 1 + 2lo)HF 

+ A 4 (67 - 2) + 37Z (33) 

= -(3~fH + F)Z. (34) 

In this analysis we discuss two cases, corresponding to Z ^ and to Z = 0. 
• In the case where 0, the first equilibrium point is given by 

F 1 = -3 7 #i (35) 
A, 



Hi = 2 — n (36) 

1 I87W - 9j 2 oj + 12 



and 



We have 



hence 



Z — ^HR. anc i p = -_3^,ip — — 3iHp 
d> a 
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in compact form 



(2A 4 + 67 A 4 + 67a; A 4 ) 
67a; — 3j 2 lo + 4 



(37) 




( J 



A 4 



V 



^\/^T87W-97 2 w+12"' 
(2A4+67A 4 +67^A 4 ) 
67W— 37 2 u+4 



• In the case where Z = (p = 0), 



w(2-2 7 ) + l\ - 



^ F 2 - ^(67 - 4)#F 



= -6 {pfoj + 1) H 2 — uj 

+ A 4 (2 7 w + 1) 
= -3 (67 - 4) H 2 - (w(4 - 37) + 3) F 2 - 3(67 - 1 + 2w)FF 
+ A 4 (6 7 -2) 

and for u; > —4/3 orw < —3/2, the equilibrium point is 
(H2, F2, Z2) : 



(38) 



2A 4 



(2lj + 3) (3w + 4) 



(a; + 1,1,0). 



(39) 



Note that for the second equilibrium point, the values of H2 and F2 are the 
same as the one in the 4-dimension case [16]. Therefore, the Z-term behaves like 
a corrective term for the 4-dimension case. In this sense we check the stability 
of this point 3 by writing 

h = H-H 2 

f = F-F 2 

z = Z-Z 2 . 

The equations (23) and (24) become 



dh 

dt 

<JL 

dt 
dz 
dt 



Ho 



The Jacobian 



( -2(3 7 w + 2) -w(27- 1) 

-9(27- 1) -(67 + 3^ + 1) 
v 



3(cj+1)7 



(3 7 + 3 7 LJ + l) J 



In the case where all the eigenvalues of the Jacobian have a non vanishing 
real part, the fixed point is called hyperbolic and the signs of this real parts 



3 Since the eigenvalues of the Jacobian of the first equilibrium point (Hi, Fi, Zi) have a 
complicated expression, its stability is not considered. 



+ 
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determine its stability. Indeed, if the real part of each eigenvalue has a negative 
sign then the equilibrium point is stable. While if the sign of the real part of 
each eigenvalue is positive, or if the sign of one of them is positive and negative 
for other, then the equilibrium point is unstable. Finally, if the real part of any 
of the eigenvalues is zero, then the equilibrium point is called nonhyperbolic 
and its stability in the neighborhood of that point cannot be determined by 
this method. 

The eigenvalues of the Jacobian at the equilibrium point (39) are: 
Ai = —(67 + 670; + !) 



2A. 



4 



(2oj + 3) (3w + 4) ' 



A 2 = -(3w + 4)< 



2A 4 



(2w + 3) (3w + 4) ' 
A 3 = -(37 + 37W + 1) 1 



2A 4 



(2w + 3) (3w + 4) 

We notice that, for 7 > and u > —4/3, the sign of this eigenvalues is 
negative and hence the equilibrium point (H2, F2, Z2) is stable. 

In this case, we should have u> > —4/3 or u < —3/2, and in the limit 
ijj — > +00 we obtain: 

H 2 ~ ~ ujF 2 (40) 

4.2 Linearized dynamical system 

To solve the dynamical system (8), (9) and (10) we linearize the solution as in 
[43]: 

H = H 2 + h(a) (41) 
F = F 2 + /(a) (42) 

Z = z{a) (43) 

where h(a), /(a) and z(a) are linearized perturbation functions to be deter- 
mined later. 

Putting (41), (42) and (43) into the field equations (8), (9), (10) and ne- 
glecting higher terms in h(a), f(a) and the product h(a)f(a) one obtains the 
following system: 



dh 
da 
d£ 
da 
1 dz_ 
\ da 



a(oj + 1) 



f -2(37^ + 2) -0,(27-1) 0±$- 2 \ 

-9 ( 27 - 1) -(67 + 30. + 1) |^ 
V -(37 + 37W + I) J 
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a 3 H 2 



I (<j(6 7 -2)+3) \ 
(2w+3) 
q (67-4) 





V 



This system becomes 



dx 
da 
dy 

da 
dz 
da 



a{u + 1) 



-3uj - 





A: 



a 3 H 2 (to + 1) 



t 

—67 — 670; — 1 


1 

(67—2^+67^—1) 
(2w+3) 








-37 — 37W — 1 



/ x 

y 

z 



(41 



with 




X + V + 3H 2 +2ujH 2 

Z 



The solutions of (45) are: 



' + 



B 

(A-2)a 2 
B> 



rA' 



-<*(*) 



(A'-2)a 2 



with 
A' 



3^ + 4. 
cj + 1 ' 

(67 + 670; + 1) 
cj + 1 
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3/e 



1 



B' = - 



H 2 (w + 1) 
3fc [67 - 2w + 670; - 1] 



#2 (w + 1) (2w + 3) 



(46) 



(47) 



and Ci, C2, C3, C[ and C 2 are integration constants. The linearized solu- 
tions (41), (42) and (43) become then: 



H = Ho- 



lt 



(u + l)(u + 3) (a \ 2 /oo\^r 
a 2 H 2 {oj + 2) (2w + 3) \a J ^ 1 \a) 



+ H K 2 



(37+37^+1) 

Ul + l 



+ 



H K 3 («L 
V a 



(67+67(^+1) 



(48) 



F 2 + 



3k 



a 2 H 2 (uj + 2) (2w + 3) V « 



a 



+ 3# ^i 



+ H «>(«> 



(37+37^+1) T ^ 

^ + 3H ^ («> 
uj \ a 



(67+67^+1) 

Ul + l 



(49) 
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, 37+37"+! 

a 1 -) — (so. 



where the subscript '0' indicates the present value. K\, K 2 and K3 are 
dimensionless integration constants. 

Letting uj — ^ 00, the linearized solutions (48), (49) and (50) are written in 
the form: 



x 3+i 



37+ ^ /n„\67+^r 



F = F 2 + 3H K 1 

+ V^P+W^P ( 52 , 
uj \ a J to \ a J 



. ., . 37+37^+1 

C 3 (i) - ,53, 

4.3 Cosmological parameters and dark energy 

In what follows, we define the individual ratios in terms of the density parameter 
p (f2j = Pi/pc) where pi run for matter, radiation, cosmological constant and 

even curvature; p c = is the critical density and Hq is the Hubble parameter 
today 

And from the standard Friedmann equations we have [44, 39]: 

Substituting the solution (51) in equation (55) and, in order to recover the 
different exponents of the equation (55), we neglect terms for which the power is 
higher than 4. Hence one gets, for each 7, in 3-brane space-time the expressions 
of the constants K\, K 2 and K% by comparing respectively the expressions of 
Qi in (55) and those of B.D cosmology in (51) for u — 1 00. 

First, let us mention that all forms of matter/energy are possible and we 
restrict ourselves to the 7-parameter of the equation of state for which the 
different exponents of the equation (55) are recovered. From (40) we have: 
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The following table summarizes the main results: 



7 


-1/3 





1/3 


1/2 


2/3 


1 


4/3 


2 


K x 




Ki 


2 v / nx 


K\ 


Um 


K x 




2 v / nx 


K 2 


V 


K 2 











V 


V 


V 


K 3 


V 


V 


V 


K z 










V 



(57) 



The character V means that all values of Ki are possible. 

In the case 7 = we have K\ + K 2 = 0, and if we take K\ = 2 ^ , then A^ 
should have the value — • And for 7 = \ and 1, we have K\ + A3 = 2 ^ A 
and if we take Ai = 2 ^ , then A3 should have the value 0. 

According to the present CMB observations and type la SNe [40, 41, 42], 
our universe seems to be spatially flat and possess a non vanishing cosmological 
constant [45]. For a flat matter dominated universe, cosmological measurements 
imply that the fraction of the contribution of the cosmological constant A 
to present energy density of the universe is JIa — 0.75 and Qm — 0.25. 

On the other hand and from the density of the microwave background pho- 
tons, pR = 4.5x10~ 34 #/cto 3 which gives Qr = 2Ah~ 2 10~ 5 where 0.4 < h < 1 
[39]. Therefore, we can safely neglect the contribution of relativistic particles 
to the total density of the universe today, which is dominated by either a non- 
relativistic particles (baryons, dark matter or massive neutrinos), a cosmological 
constant or an exotic form of matter /energy. 

An interesting consequence of these considerations is that one can write the 
Friedmann equation (55) today as: 



i = n A + n M + n k (58) 

In what follows, we discuss all possible form of matter /energy, so that we 
recover the different exponents of the equation (55). 



4.3.1 Flat universe 

According to equation (58), the line f^A = 1 — VLm corresponds to a flat universe 

(Qk = 0), and separates the open universe from the closed one. 

Except the cases where 7 = 0, \ and 1, the table (57) shows that K\ = . 

If K\ = 27j7§p the universe becomes flat. If K\ < , we obtain an open 
universe, otherwise the universe is close. 

With the numerical value, K\ ~ 0.144 34, we conclude that our universe is 
flat and the theory is in agreement with the observation data. 



4.3.2 Accelerating universe 

Consider now the deceleration parameter [3, 44, 39] 

q = tt k + ^tt M - tt A . (59) 

Using the present results obtained on density parameters we neglect Q k and 
one can parameterize the matter /energy content of the universe with just two 



14 



components: the matter, characterized by and the vacuum energy charac- 
terized by i.e., qo = \£Im — ^A- 

A uniform expansion (go = 0) corresponds to the line = separating 
the accelerating from the decelerating universe and K\ verify: 



If K\ < \/TT\, the universe is in an accelerating phase while K\ > \ZTF\ corre- 

sponds to a decelerating phase of the universe. 

Consequently, the (Qm,^a) plane shows that we live in an accelerating 
flat universe, since numerical calculations show that K\ < y/Tl\, which is in 
accordance with the experimental data of la SNe [40]. 

5 Conclusions 

In this work we have examined the behavior of 3-brane of Brans-Dicke cosmol- 
ogy which differs from other Brans-Dicke cosmology by the fact that the 5D 
approach affects the ordinary matter (by the square of energy density) but not 
the Brans-Dicke field. This approach clearly shows how to describe the early 
and late time behavior of the universe first by means of the scale factor (section 
3) and second by means of the cosmological parameters (section 4) for large 
value of the a;-parameter. Let us notice also that the present work is a general- 
ization of our previous work [16] in which we did not consider the intermediate 
case. 

Furthermore, this approach gives two possibilities to describe the universe. 
The first one, consider that the universe underwent a bounce state and hence 
avoided the big bang singularity. The second one in which the universe will 
undergo a turnaround state and therefore avoiding, probably, the big rip sin- 
gularity. In the other cases where no bounce nor turnaround state are present, 
the universe began to expand from the big bang singularity in the past, or/and 
will meet the big rip singularity and its dramatic consequences. The two pos- 
sibilities are consistent with the fact that, today, the universe undergoes an 
accelerating period. However, we can opt for the case in which the universe 
underwent a bounce state in the past and will undergoe a turnaround state in 
the future in order to avoid the dramatic consequences of the big rip and to 
have the possibility of crossing the phantom divide line. 

Finally, we conclude that the assumption of 3-brane behavior of Brans-Dicke 
cosmology gives an interesting results and enables us to explore this approach 
in more detail in future investigations. 
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